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ABSTRACT 

We show that if FU Orionis objects harbour hot Jupiters embedded in their discs, 
the resulting non-axisymmetric dissipation profile in the disc would be manifest as 
time-dependent distortions in the absorption line profiles of these objects. In order 
to affect the infrared line profiles, planets must lie within ~ 0.5 au of the central 
star, whereas only planets within ~ 0.1 au would influence the optical line profiles. 
The timescale for modulation of the line profiles is relatively short (months) in each 
case, so that the effect could not have been discovered from published spectra (which 
combine data taken in different observing seasons). The detection of hot Jupiters in 
FU Orionis objects would be in line with the expectations of tidal migration theories 
(which predict a high incidence of close planets around young stars) and would also 
lend support to models which link the triggering of rapid rise FU Orionis events to 
the existence of a close massive planet. 

Key words: accretion, accretion discs — line: profiles — planets and satellites: 
formation — planetary systems: protoplanetary disks — stars: pre-main sequence — 
stars: individual (FU Orionis) 



1 INTRODUCTION 

Radial velocity surveys show that the frequency of massive 
planetary companions to late F and G dwarf stars on the 
main sequence is at least 6 percent. Around 1 percent of 
surveyed stars possess 'hot Jupiters', that is, massive plan- 
ets within a tenth of an au of their parent star (Butler et al. 
2000; Marcy & Butler 2000). Planetary migration theories 
(Lin & Papaloizou 1986; Lin, Bodenheimer & Richardson 
1996) suggest that the abundance of such planets could be 
considerably higher in pre-main sequence stars, but even 
the main sequence numbers imply that a significant frac- 
tion - though probably not a majority - of pre-main se- 
quence (T Tauri) stars harbour planetary companions at rel- 
atively small orbital radii. Directly detecting planets around 
T Tauri stars is, however, problematic. The rapid rotation 
of young stars means that searches for planetary transits 
have to contend with additional variability due to rotational 
modulation on a timescale of a few days (Aigrain, private 
communication), whilst the strong veiling of stellar absorp- 
tion features in Classical T Tauri stars have made it difficult 
to detect even stellar companions by radial velocity methods 
(Mathieu 1994). 



In this paper, we point out that the subset of accret- 
ing pre-main sequence stars known as FU Orionis objects 
present a novel opportunity for searching for hot Jupiters. 
FU Orionis objects are believed to be episodes of enhanced 
accretion on to young stars, during which the emission 
from the central star is completely swamped by the ac- 
cretion luminosity (corresponding to accretion rates of ~ 
10 -5 — lO _4 M0yr _1 ; see Hartmann, Kenyon & Hartigan 
1993; Kenyon 1995; Hartmann & Kenyon 1996 for reviews). 
During the outburst, both the broad band spectral energy 
distribution and the broad double peaked absorption lines of 
FU Orionis objects are well fit by the predictions of a steady 
state Keplerian accretion disc (Hartmann & Kenyon 1985, 
1987; Kenyon, Hartmann & Hewett 1988; Popham et al. 
1996). In Section 2, we show that the presence of a planet in 
the inner disc of an FU Orionis object would generate a non- 
axisymmetric pattern of disc emission, due to a high rate of 
accretion onto the planet from the disc. In consequence, the 
resulting line profiles are likely to be asymmetric and time 
variable, with the extra absorption moving between the red 
and blue wing of the line as the planet orbits. We present 
simple models of the predicted line profiles in this case. In 
Section 3 we discuss which lines provide the best diagnostics 
of hot Jupiters, and discuss the implications of a detection 
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for planetary migration theory and the origin of FU Orionis 
outbursts. 



2 NON-AXISYMMETRIC DISSIPATION AND 
SPECTRAL LINE PREDICTIONS 

The interaction of a planet embedded within a viscous ac- 
cretion disc has several facets. First, if the planet is suffi- 
ciently massive, gravitational torques from the planet will 
clear a gap (an annular region in which the surface den- 
sity is greatly reduced as compared to the disc) in the disc 
(Goldreich & Tremaine 1980; Papaloizou & Lin 1984). Pa- 
rameterizing the viscosity v using the Shakura & Sunyaev 
(1973) Q-prescription, v = ac^/Q, where c s is the sound 
speed and ft the angular velocity, gap opening occurs for a 
mass ratio q = M p /M* given by (Takeuchi, Miyama & Lin 
1996), 



5crit ; 



1/2 



(1) 



where h/r is the relative scale height of the disc. In the case 
of FU Orionis outbursts, detailed time-dependent modeling 
of the disc suggests that h/r ~ 0.2 (Bell & Lin 1994), while 
numerical simulations of angular momentum transport in 
discs typically find a ~ 10" 2 (e.g. Stone et al. 1996). The 
stellar mass is probably comparable to a Solar mass or some- 
what smaller - for FU Orionis itself Hartmann & Kenyon 
(1987) infer M* ~ 0.6 Mq. Using these parameters, the crit- 
ical value for q is around 4 x 10~ 3 . This estimate suggests 
that because of the increased disc thickness during FU Ori- 
onis events, Jupiter mass planets around Solar mass stars 
will no longer be able to open a clean gap, and will instead 
be embedded within disc gas (Clarke & Syer 1996). 

Second, irrespective of the presence or absence of a gap, 
there can be ongoing accretion from the disc onto the planet 
(Artymowicz & Lubow 1996). The planetary accretion rate 
can be expressed as a fraction of the disc accretion rate, 



M p = e(q)M dil 



(2) 



where Mdi sc is the accretion rate in the disc at a radius sub- 
stantially greater than that of the planet. Simulations by 
Lubow, Seibert & Artymowicz (1999), D'Angelo, Henning 
& Kley (2002), and Bate et al. (2002) show that the effi- 
ciency factor e ~ 1 for a planet that is marginally able to 
open a gap. For larger masses, there is an approximately 
exponential suppression of the accretion efficiency (Lubow, 
Seibert & Artymowicz 1999). For lower mass planets e also 
decreases, with a scaling variously measured as M^ 3 from 
2D simulations (D'Angelo, Henning & Kley 2002) and as 
M p in 3D calculations (Bate et al. 2002). In either event, 
Jupiter mass planets, embedded within the inner disc of FU 
Orionis objects, are expected to accrete mass at a significant 
fraction of the (very high) disc accretion rate. 

Finally, the interaction of massive planets with the disc 
may induce eccentricity in both the disc material and the 
planetary orbit. Whether this occurs for Jupiter mass plan- 
ets is currently uncertain (Papaloizou, Nelson & Masset 
2001; Goldreich & Sari 2003; Ogilvie & Lubow 2003), and 
is probably dependent upon details of the angular momen- 
tum transport process in protoplanetary discs (e.g. Ogilvie 
2001). We will ignore this possibility henceforth. 




Figure 1. Surface density from an isothermal, two dimensional 
hydrodynamic simulation of the planet-disc interaction, with pa- 
rameters appropriate for an FU Orionis event. We modeled a 
planet of mass 3 Mj, on a circular orbit in a disc with Cs/v^ = 
0.15 at the orbital radius of the planet. Angular momentum trans- 
port was approximated using a constant kinematic viscosity cor- 
responding to an a of 4 X 10~ 3 at the radius of the planet. The 
surface density is plotted on a linear scale, with S = 1 (in arbi- 
trary units) corresponding to the initial value of the disc surface 
density at large radius. To emphasize the structure in the disc, 
high values of S close to the planet are suppressed in the image. 
Scales on the x and y axes are in units of the planetary orbital ra- 
dius. As discussed in the text, we consider orbital radii r ~ 0.1 au, 
though this physical scale is not an input for the simulation. 



Planet mass 3 Mj 

Stellar mass 1 Mq 

Implied disc thickness (at planet radius) (h/r) ~ 0.2 

Effective a (at planet radius) a = 4 X 10~ 3 

Table 1. Main parameters of the system shown in Figure 1. Note 
that the simulation can be rescaled to represent other systems 
with the same mass ratio and h/r. Likewise, the planetary orbital 
radius and disc accretion rate used are arbitrary. 



To verify some of these estimates, we have simulated the 
interaction between the planet and the disc using the ZEUS 
hydrodynamics code (Stone & Norman 1992). The most im- 
portant parameters for the disc and planet are summarized 
in Table 1. The simulation was run at a numerical resolution 
of 240 2 grid points in cylindrical co-ordinates, with outflow 
boundary conditions at the radial edges of the grid. After 
120 orbits of the embedded planet, by which time initial 
transients have died away, the surface density is as shown in 
Figure 1. As expected based on the arguments given above, 
the planet (here of mass 3 Mj, where Mj is the mass of 
Jupiter) is unable to clear a well-defined gap in a disc with 
relative scale height h/r appropriate to FU Orionis systems, 
though there is a clear drop in the surface density at the ra- 
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dius where a gap would form in a cooler disc. A stream from 
the outer disc carries approximately half of the outer disc 
accretion rate, forming an (unresolved) circumplanetary disc 
and extending toward the inner boundary of the simulation. 

Apart from gas in the (small) circumplanetary disc, the 
velocity perturbations to the disc flow are negligible. We 
have therefore modeled the influence of the planet upon the 
line profiles by assuming that the dominant effect is due 
to the enhanced dissipation caused by accretion onto the 
planet. The magnitude of the extra dissipation is, 



AL 



GM p M diB 



(3) 



where R p is the planetary radius. Since the disc is extremely 
optically thick, this extra luminosity will be radiated over a 
hotspot with a radius ~ h ~ 0.2r. The ratio of the effective 
temperature of the hotspot T e to the effective temperature 
T e of the disc is then, 



2? 



L + AL 



(4) 



where L is the luminosity of the hotspot area in the absence 
of the enhancement due to the planet. This unperturbed 
luminosity is given by the product of the hotspot area (2Tvh 2 , 
allowing for both sides of the disc) and the radiative flux 
generated by a steady-state accretion disc (e.g. Frank, King 
& Raine 1992). This product yields, 



L = 



3h GM*M dis 
4r 3 



(5) 



From equations @ and JSJ, it follows that, 
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where Rj is the radius of Jupiter. As noted earlier, a disc 
aspect ratio of h/r « 0.2 is consistent with detailed modeling 
of FU Orionis events as thermal disc instabilities (Bell & Lin 
1994). The planetary radius is more uncertain. In isolation, 
giant planets shrink rapidly - a 3Mj planet contracts below 
27? j within a few Myr (Burrows et al. 1997). A planet which 
is still actively accreting will be larger, so for this estimate 
we take a value of 5Rj. 

For the parameters given in equation we estimate 
that the dissipation rate per unit area is about a factor of 
two greater in the vicinity of the planet than at other phases. 
For optically thick emission, this corresponds to a 20 per- 
cent increase in the effective temperature of the disc in the 
region around the planet. In the case of absorption lines 
whose equivalent width varies only weakly with tempera- 
ture (see below) this extra dissipation produces additional 
line absorption in the velocity channels centred on the in- 
stantaneous line of sight velocity of the planet. We note that 
the resulting local enhancement of the disc temperature does 
not depend upon the accretion rate through the disc, and 
(for fixed planet radius) becomes more significant at larger 
disc radii. 

In Figures 2 and 3 we illustrate the distortion of the 
line profiles resulting from such a non-axisymmetric surface 
brightness distribution. In each case, we have plotted the 
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Figure 2. Rotational broadening function for absorption lines in 
the spectral region centred on 2.2 microns (see text for details). 
The normalisation is arbitrary and half the line is plotted (the 
other half of the line would be the mirror image of the unper- 
turbed profile, shown by the solid histogram). The dashed his- 
tograms correspond to models in which a 'planet' is placed at the 
radii indicated. 
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Figure 3. As for Figure 2, except that we show the rotational 
broadening function for the spectral region centred on 6170 A. 
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instantaneous line profiles when the planet's velocity is most 
closely aligned to the observer's line of sight. Half the line 
is plotted, the other half of the line being the mirror image 
of the unperturbed profile shown by the solid histogram. 
As the planet orbits the central star, the feature in the line 
profile would migrate back and forth between the red and 
blue wing of the line on the planet's orbital period. 

The quantity plotted is the rotational broadening func- 
tion (with arbitrary normalisation) i.e. it is the function 
with which the intrinsic spectrum would be convolved in the 
relevant region of the spectrum. In practice, the rotational 
broadening function in observed spectra of FU Orionis ob- 
jects is extracted using cross-correlation techniques, which 
have the twin advantage of providing a weighted average 
over many lines and at the same time being mainly sensitive 
to line shape, rather than line depths (see Hartmann et al. 
1986). In the construction of Figures 2 and 3 we adopt a sim- 
plified model in which the absorption per unit area at each 
point in the disc is simply proportional to the black body 
emission flux at that point, (i.e. we neglect variations in the 
line equivalent width with either temperature or gravity over 
temperature ranges appropriate to the lines in question). A 
more sophisticated treatment would involve the summation 
of lines from appropriate stellar models, although it is un- 
clear whether this is in fact more correct, given that in discs 
there is likely to be significant deposition of mechanical en- 
ergy in the atmosphere, which makes comparison with stel- 
lar spectra problematical. In any case, Kenyon et al. (1988) 
noted that in the infrared, the saturated nature of the CO 
lines makes the equivalent widths rather insensitive to the 
effective temperature and gravity and a similar argument 
can be made for the principal lines in the 6170 A window 
(see discussion in Hartmann and Kenyon 1987). 

The two spectral regions illustrated in Figures 2 and 3 
correspond to bands centred on 2.2/wn and 6170 A, and, fol- 
lowing Hartmann and Kenyon 1987 and Kenyon et al 1988, 
we have considered contributions to the line profiles in each 
spectral region from the temperature ranges 900 — 5300 K 
and 3300 — 6600 K respectively. For the temperature pro- 
file of the unperturbed disc, we have adopted the model of 
Kenyon et al. (1988) which provided the best fit to the spec- 
tral energy distribution and line profiles of V1057 Cygni: 
the disc temperature follows the classic r~ 3//4 scaling with 
radius, with maximum temperature of ~ 9000K at the disc's 
inner edge (0.016 au). Such a disc temperature distribution 
corresponds to an accretion rate of ~ lO~ 4 Af0yr -1 onto a 
solar mass star. Note that material at the inner edge of the 
disc is too hot to contribute to the absorption in either of 
the wavebands shown here and so there is likewise no con- 
tribution from any (presumably hotter) boundary layer. 

In the absence of the planet, the absorption lines exhibit 
the double peaked structure familiar from the emission line 
profiles of cataclysmic variables. We have here scaled the ve- 
locities so as to roughly reproduce the line widths in V1057 
Cyg, requiring the projected Keplerian velocity at the stellar 
surface to be ~ 40 km s _1 (Kenyon et al. 1988). The max- 
imum velocity corresponds to the projected velocity at the 
inner edge of the line producing region, whereas the velocity 
at the line peak reflects the projected velocity at the location 
where the line emissivity starts to fall (close to the Wien cut- 
off of the relevant transition). At lower velocities, the line is 
produced by a mixture of high velocity material projected 



Figure 4. Grey scale representation of the evolution of the ro- 
tational broadening function for the spectral region centred on 
2.2/im. The horizontal axis is velocity over the range —8 km/s (left 
hand margin) to +46 km/s (right hand margin). Each horizon- 
tal slice is constructed by subtracting the mean (orbit-avcraged) 
value of the rotational broadening function from its instantaneous 
value in that velocity channel. (The black denotes zero difference, 
whilst the whitest pixel corresponds to a residual equal to about 
10 percent of the orbit averaged value of the function in that ve- 
locity channel). The spectrum has been degraded to a resolution 
of 4 km/s and the vertical (time) axis corresponds to 3/8 of the 
orbital period of the planet. 



at low velocity and by low velocity material absorbing in the 
Wien region of the spectrum. For the parameters employed 
here, the regions of the unperturbed disc that contribute 
to the line absorption lie between 0.025 — 0.07 au (for the 
optical lines) and 0.04 — 0.4 au (for the infrared lines). 

We model the presence of the planet as a 20 percent 
boost in the effective temperature distributed in a region of 
azimuthal extent 0.17T and radial half-width 0.1r p centred 
on the planet at radius r p , these choices being motivated 
by the preceding discussion. The effect of the planet is thus 
manifest in the absorption line profiles as additional absorp- 
tion at the projected Keplerian velocity of the planet. It is 
evident from Figures 2 and 3 that the planet produces a dis- 
tortion to the line profile whenever its associated region of 
enhanced dissipation overlaps the region of the disc produc- 
ing the relevant line. The feature becomes less conspicuous, 
however, at lower velocities (i.e. within the peak of the line 
profile), since the contribution from the region disturbed by 
the planet is then diluted by the contribution from higher 
velocity material projected into this spectral region. 

We note that the distortion of the line profile is quite 
marked. For example, for the planet located at 0.25 au, the 
extra emission from the vicinity of the planet represents 
around 22% of the continuum at 2.2/im, produced in a region 
of the disc from which there is significant line absorption at 
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this wavelength and in a velocity channel around 10 km/s. 
Thus the line absorption in this channel is boosted by this 
factor. The extra emission from around the planet however 
only boosts the overall continuum level at 2.2/xm by 1.5% so 
that the line to continuum ratio in this channel is boosted by 
nearly 20%, as shown in Figure 2. We note that the small 
contribution of the planet to the over-all continnum level 
means that such a planet would not produce any detectable 
effect on the broad band spectral energy distribution of the 
system. 

The features generated by the planet in Figures 2 and 3 
are narrow, with velocity widths - comparable to the Kep- 
lerian shear across regions of the disc strongly perturbed by 
the planet - of only a few km/s. As such, they would be hard 
to resolve in individual spectra, even given the high spectral 
resolution of a new generation of infrared spectrometers such 
as NIRSPEC and Phoenix. On the other hand, line mon- 
itoring should reveal these features through their regular 
motion across the line profile. Figure 4 is a grey scale repre- 
sentation of the spectral variations produced at 2.2^im by a 
planet orbiting at 0.125 au, viewed at a spectral resolution 
of R=75,000. Each horizontal slice is a spectrum from which 
the mean line profile (over an orbit of the planet) has been 
subtracted. The planet's motion is along the line of sight at 
the bottom of the plot and the feature moves through the 
line centre (left hand margin) as the planet motion passes 
through the tangent plane. The characteristic trail produced 
by the planet is evident despite the inconspicuousness of the 
planet's feature in any individual spectrum. 



3 THE DETECTABILITY OF HOT JUPITERS 
IN FU ORIONIS SYSTEMS 

The above crude model suggests that the sort of non- 
axisymmetric dissipation profile that one would associate 
with an embedded hot Jupiter would have a conspicuous ef- 
fect on the line profiles of FU Orionis systems. This immedi- 
ately raises the question of whether such a distortion should 
already have been detected in FU Orionis objects analysed 
to date if they indeed harboured such planets. Two factors 
however militate against the discovery of such features in 
existing spectra: firstly the limited spectral resolution (par- 
ticularly in the infrared, where the spectral resolution of 
published spectra is comparable with the line width) and 
secondly, the fact that in order to enhance the signal to 
noise of the lines analysed, spectra have been co-added from 
a number of observing runs. We have seen above that the 
regions of the disc in which embedded planets would be de- 
tectable in this way are close in (within 0.5 au) and hence 
the associated orbital timescale is short (months). The effect 
of a feature crossing the line profile from the red to the blue 
wing would thus be washed out by combining spectra from 
different observing seasons. 

In order to detect an embedded planet through its time- 
dependent distortion of the line profile, it is of course nec- 
essary that the planet be located in the region of the disc 
producing the line. This restricts the technique to the de- 
tection of rather close planets (within ~ 0.5 au for an effect 
to be manifest in the infrared lines and within ~ 0.1 au for 
a manifestation in optical line spectra). Such orbital radii 
are well within the range of separations of hot Jupiters (to 



date 37 and 18 exoplanets respectively have been discovered 
in each of these separation ranges 1 ). Since the infrared lines 
are produced in cooler regions of the disc, it follows that 
the infrared lines provide a good diagnostic of planets over 
a larger range of radii. On the other hand, since the optical 
lines originate in a more spatially restricted region of the 
disc, it is evident that if a, planet is present at an appropri- 
ate radius, it will be more conspicuous in the optical line 
profiles, since the distortion produced by the planet is less 
diluted by absorption by other annuli. 

What, then, is the likelihood of there being an appro- 
priate planet in a FU Orionis system? If the incidence of 
planets in FU Orionis systems follows the statistics of ex- 
oplanets around main sequence stars, then we expect that 
between 2 and 3 percent of systems will harbour a planet 
within 0.5 au. There would then be little chance of there 
being an exoplanet in the « 16 systems classified as FU 
Orionis systems to date (Sandell & Weintraub 2001). There 
are, however, strong theoretical arguments that close planets 
should have been commoner during the pre-main sequence 
stage, and that the incidence of such systems on the main se- 
quence represents only the (possibly small) fraction that sur- 
vived migration without being swallowed by the star (Lin, 
Bodenheimer & Richardson 1996; Trilling et al. 1998; Ar- 
mitage et al. 2002; Trilling, Lunine & Benz 2002). Indeed, 
given the extreme youth of most FU Orionis systems, it is 
probable that all close planets present in FU Orionis discs 
are destined to be swallowed, in which case main sequence 
statistics provide no constraint whatsoever on the possible 
population. Although it is impossible to assess by what fac- 
tor the incidence of hot Jupiters is increased in young stars, 
the sign of the effect should motivate the search advocated 
here. 

Finally, we remark that FU Orionis objects may be 
prime targets for detecting hot Jupiters if, as suggested by 
Clarke and Syer (1996), there is a causal link between the 
two phenomena. Clarke and Syer (1996) argued that the 
most plausible way of triggering a rapid rise outburst (i.e. 
one in which the hundredfold increase in luminosity occurs 
on a timescale of a few months), is if the disc prior to out- 
burst is 'dammed up' behind such a planet. In the Clarke and 
Syer model, the onset of the thermal ionisation instability 
upstream of the planet causes the dam to burst, enveloping 
the planet and causing a rapid rise in accretion luminosity. 
It is notable that the planetary characteristics advocated by 
Clarke and Syer (a 10 Jupiter mass planet located at ~ 0.1 
au), are just such that they would be detectable through 
the absorption line diagnostics discussed here. V1057 Cygni 
and FU Orionis (both of which are of the rapid rise variety 
and which are well studied spectroscopically) are thus both 
prime candidates for this sort of investigation. 
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